Measurement of the Flux of Ultrahigh Energy Cosmic Rays from Monocular 
Observations by the High Resolution Fly's Eye Experiment. 
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We have measured the cosmic ray spectrum above 10 17 ' 2 eV using the two air fluorescence detec- 
tors of the High Resolution Fly's Eye observatory operating in monocular mode. We describe the 
detector, photo-tube and atmospheric calibrations, as well as the analysis techniques for the two 
detectors. We fit the spectrum to a model consisting of galactic and extra-galactic sources. 



The highest energy cosmic rays detected so far, of ener- 
gies up to and exceeding 10 20 eV, are very interesting in 
that they shed light on two important questions: the na- 
ture of their origin in astrophysical or other sources and 
their propagation to us through the Cosmic Microwave 
Background Radiation (CMBR) . The production of pions 
from interactions of CMBR photons and Ultra High En- 
ergy Cosmic Rays (UHECR) is an important energy loss 
mechanism above ~ 10 19 ' 8 eV, and produces the Greisen- 
Zatsepin-K 'uzmin (GZK) effect 0,0; e+e" production 
in the same collisions is a weaker energy-loss mechanism 
above a threshold of 10 17 8 eV. We report here the flux 
of UHECR from 10 17 2 eV to over 10 20 eV, measured 
in monocular mode, with the High Resolution Fly's Eye 
(HiRcs) detectors. 

The HiRes observatory consists of two air-fluorescence 
detector sites separated by 12.6 km and located at the 
U.S. Army Dugway Proving Ground in Utah. Cosmic 
rays interacting in the upper atmosphere initiate parti- 
cle cascades known as extensive air-showers. Passage of 
charged particles excites nitrogen molecules causing emis- 
sion of (mostly) ultraviolet light. The fluorescence yield 
has been previously measured by Kakimoto et al. 0] , and 
more recently by Nagano et al. J4| For this analysis, we 
used the fluorescence spectrum compiled by Bunner |5j 
and normalized it to the yield of Kakimoto. By mea- 
suring the longitudinal development of the fluorescence 
signal, one can infer the arrival direction, energy and av- 



erage composition of the primary cosmic ray. HiRes was 
designed to measure the fluorescence light stereoscopi- 
cally. However our two detectors trigger and reconstruct 
events independently. In this "monocular" mode our cur- 
rent data have significantly better statistical power and 
cover a much wider energy range than our stereo sample. 

The two HiRes detector sites, referred to as HiRes- 
I and HiRes-II, are operated on clear, moon-less nights. 
Over a typical year, each detector accumulates up to 1000 
hours of observation. The HiRcs-I site has been in op- 
eration since June of 1997 @|. It consists of 21 detector 
units, each equipped with a 5 m 2 spherical mirror and 
256 photo-tube pixels at its focal plane. Each photo- 
tube covers a 1° cone of sky. These 21 mirrors cover 
elevation angles between 3° and 17°. The HiRes-I elec- 
tronics perform sample-and-hold integration in a 5.6 [is 
window, which is long enough to contain signals from all 
reconstructible events. The HiRes-II site was completed 
in late 1999 and began observations that year. This site 
uses the same type of mirrors and photo-tubes as HiRes- 
I, but contains 42 mirrors, in two rings, covering elevation 
angles from 3° to 31°. HiRes-II uses an FADC data acqui- 
sition system operating at 10 MHz [7|. Both the HiRes-I 
and -II sites provide 2-7T azimuthal angle coverage. 

To determine the correct shower energies, the air flu- 
orescence technique requires accurate measurement and 
monitoring of photo-tube gains. Two methods of calibra- 
tion are used. Pulses from a YAG laser are distributed to 
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mirrors via optical fibers. They provide a nightly relative 
calibration. A stable, standard light source is used for a 
more precise monthly absolute calibration. Overall, the 
relative photo-tube gains were stable to within 3.5% and 
the absolute gains were known to ±10% 0- 

A second variable in the energy measurement is at- 
mospheric clarity. Light from air showers is attenuated 
by: (a) molecular (Rayleigh), and (b) aerosol scatter- 
ing. The former is approximately constant, subject only 
to small variations in the atmospheric overburden. The 
aerosol concentration varies with time. HiRes measures 
the aerosol content by observing scattered light from two 
stecrablc laser systems. The laser observed by HiRes-I 
has been in operation since 1999. The vertical aerosol 
optical depth from ground to 3.5 km altitude, ta, is 
measured each hour (the vertical transmission through 
the aerosol is T = e~ TA ). Over two years, these mea- 
surements yielded an average ta at 355 nm of 0.04 0. 
The RMS of the distribution is 0.02, and the systematic 
uncertainty in the mean is no larger than this. The av- 
erage aerosol ground-level horizontal extinction length, 
Ah, was determined to be 25 km. Because about half of 
our data were collected before the steerable lasers were 
in operation we used these averages in our analysis and 
simulation. 

Between June 1997 and February 2003, the HiRes-I de- 
tector operated for approximately 3600 hours. From this, 
2820 hours of good weather data were analyzed. We se- 
lected 5.5 million downward, track-like events. For each 
of these, a shower-detector plane was determined from 
the pattern of photo-tube hits. We excluded events con- 
taining an average number of photo-electrons per photo- 
tube of less than 25, where fluctuations in signals are too 
great to permit reliable reconstruction. We also cut out 
tracks with angular speed in excess of 3.33°//^s; for these 
events (typically within 5 km of HiRcs-I) the shower max- 
ima appear above the field of view. We selected 12,709 
events for reconstruction. 

Determination of the shower geometry is possible in 
monocular mode. The impact parameter, R p , and the 
angle of the shower in the plane containing the shower 
and the detector, ip, are found by fitting the photo-tube 
trigger times to the angles at which they view the shower. 
However, HiRes-I monocular events are too short in an- 
gular spread for reliable pure-timing fit. For this anal- 
ysis, the expected form of the shower development was 
used to constrain the time fit to yield realistic geome- 
tries. The shower profile is assumed to be described 
by the Gaisser-Hillas parameterization ^lj, which has 
been found to be in good agreement with previous HiRes 
measurements and with CORSIKA/QGSJET simu- 
lations [T3L H4I Il5| . This technique is called the Profile- 
Constrained Fit (PCF). We allowed the shower maxi- 
mum, x m to vary in 35 g/cm 2 steps between 680 and 
900 g/cm 2 , matching the expected range of x m for pro- 
ton to iron primaries. After reconstruction, we require a 
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FIG. 1: Comparison of HiRes-I simulated (histogram) and ob- 
served (points) R p distributions at (a) 10 18 ' 5 , (b) 10 190 , and 
(c) 10 19 ' 5 eV. The Monte Carlo distributions are normalized 
to the number of data events. 



minimum track arc length of 8.0° and a maximum depth 
for the highest elevation hit of 1000 g/cm 2 . Significant 
contamination from the forward-beamed direct Cerenkov 
light degrades the reliability of the PCF. Therefore, we 
rejected tracks with ip > 120° and those with two or more 
angular bins of the shower with >25% Cerenkov light. A 
total of 6,920 events remained. 

Monte Carlo studies were performed to assess the re- 
liability of the PCF method. The simulated events were 
subjected to the same selection criteria and cuts imposed 
on the data. An RMS energy resolution of better than 
20% was seen above 10 19 5 eV. However, the resolution 
degrades at lower energies to about 25% at 10 18 5 eV. 
These Monte Carlo results were cross-checked by exam- 
ination of a smaller set of stereo events where the ge- 
ometry is more precisely known. Comparing the energies 
reconstructed using monocular and stereo geometries, we 
obtained resolutions similar to those seen in simulation. 

The simulation is also used to calculate the aperture. 
To verify the reliability of this calculation, we compared 
Monte Carlo distributions of many geometrical and phys- 
ical variables to the actual data, and consistently found 
good agreement. The Monte Carlo predictions for the 
zenith angle and impact parameter (R p ), in particular, 
are sensitive to the detector operating parameters. We 
use input parameters representative of actual running 
conditions, and again see good agreement between data 
and simulation. For example, we show the comparison of 
R p distribution at three energies in Figure^ 

The analysis of HiRes-II monocular data was similar 
to that for HiRes-I. The data sample was collected dur- 
ing 142 hours of good weather between Dec. 1999 and 
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FIG. 2: A comparison of the number of photo-electrons 
per degree of track seen in HiRes-II monocular events (data 
points) and in simulation (histogram). 



May 2000. This period represents the first stable run- 
ning of the HiRes-II detector. At the end of this period, 
a considerable change was made in the trigger, so that 
subsequently collected data will be analyzed separately. 
With the greater elevation coverage at HiRes-II, it was 
feasible to reconstruct the shower geometry from timing 
alone (the PCF is unnecessary). Therefore we were able 
to loosen some cuts for the HiRes-II fits. At this stage 
104,048 downward-going events remained. 

With the geometry of the shower known, we fit the 
measured shower profile to the Gaisser-Hillas parameter- 
ization |llj |. The events were required to have a good 
fit to the Gaisser-Hillas function, to have a track length 
greater than 10° for upper ring or multi-mirror events, a 
track length greater than 7° for lower ring events, an an- 
gular speed less than ll°/jus (the larger cut for HiRes-II 
reflects its extended elevation coverage), a zenith angle 
less than 60°, and a shower maximum visible in our de- 
tector. There was also a cut on the size of the Cerenkov 
light subtraction at < 60% of signal. Again, the same se- 
lections and cuts were applied to both simulated and real 
events. There were 781 events left after these cuts. These 
simulations also gave excellent reproduction of the data, 
as seen, for example, in the comparison of the number of 
photo-electrons per degree of track in Fig. [21 

For both HiRcs-I and HiRes-II events, the photo- 
electron count was converted to a shower size at each 
atmospheric depth, using the known geometry of the 
shower, and corrected for atmospheric attenuation. We 
integrated the resulting function over x and then multi- 
plied by the average energy loss per particle to give the 
visible shower energy. A correction for energy carried 
off by non-observable particles to give the total shower 
energy (~ 10%) was then applied. 

The monocular reconstruction apertures are shown in 
Fig. 03 Both approach 10 4 km 2 sr above 10 20 eV. We re- 
strict our result for HiRes-I to energies > 10 18 ' 5 eV; below 
this the PCF technique is unstable. Due to longer tracks 
and additional timing information, the RMS energy res- 
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FIG. 4: Combined HiRes monocular spectrum. The squares 
and circles represent the HiRes-I and II differential flux J(E), 
multiplied by _E 3 . The error bars are statistical only, and the 
systematic uncertainties are indicated by the shaded region. 
The line is a fit to the data of a model, described in the text, of 
galactic and extra-galactic cosmic ray sources. The AGASA 
spectrum [l6j| is shown by triangles for comparison. 



olution for HiRes-II remains better than 30% down to 
10 17 ' 2 eV. However, the HiRes-II data become statisti- 
cally depleted above 10 19 eV. 

We calculated the cosmic ray flux for HiRes-I above 
10 18 - 5 eV, and for HiRes-II above 10 172 eV. This com- 
bined spectrum is shown in Fig. 0] where the flux J(E) 
has been multiplied by E 3 . The error bars represent the 
68% confidence interval for the Poisson fluctuations in 
the number of events. The HiRcs-I flux is the result of 
two independent analyses [lT], H3 • 

The largest systematic uncertainties in the energy 
scale are the absolute calibration of the photo-tubes 
(±10%) @, the fluorescence yield (±10%) 0, fl, and 
the correction for unobserved energy in the shower 
(±5%) [HG3. Excluding atmospheric effects, the en- 
ergy scale uncertainty is ±15%. This translates to a sys- 
tematic uncertainty in the flux, J(E), of ±27%. 
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We estimate the atmospheric contribution to the en- 
ergy error by repeating the event reconstruction with ta 
varied by ±1 RMS value, from 0.04 to 0.06 and 0.02. 
While a ±0.02 change in ta represents a ±50% change 
in aerosol concentration, that contribution to the attenu- 
ation at these levels is small, and a 0.02 change modifies 
the transmission by only ~ 10% at 25 km from the detec- 
tor. We found the reconstructed geometries of HiRes-I 
events above 10 18 5 eV to be insensitive to changes in as- 
sumed ta, and we saw a maximum change in the energy 
of ±13% at 10 20 eV, decreasing to ±6% at 10 185 eV. The 
geometries of the HiRes-II events do not depend on ta at 
all, and the change in energy scale for these are typically 
6% or less. Taking the HiRes-I average energy shift, 9%, 
the overall systematic uncertainty in energy scale, includ- 
ing atmospheric effects, then becomes ±17%. 

We also recalculated the aperture and the flux with ta 
changed by ±0.02. From these we obtained an average 
atmospheric uncertainty in J(E) of ±15%, and an uncer- 
tainty in the flux of ±31%. The overall systematic un- 
certainties in the flux, including the energy-dependent at- 
mospheric contribution (with a slope in \A.J(E)\ of about 
5% per decade of energy) are indicated by the shaded re- 
gion in Figure 0| The relative calibration uncertainty 
between the two detector sites is less than 10%. 

Our data contain two events at or above 10 20 eV, mea- 
sured at 1.0 and 1.5 x 10 20 eV. Assuming a purely molecu- 
lar atmosphere (ta — 0.0), we obtain lower energy limits 
of 0.9 and 1.2 x 10 20 eV. In the energy range where both 
detectors' data have good statistical power, the results 
agree with each other very well. However, our flux values 
are on average 13% lower than the stereo spectrum re- 
ported by Fly's Eye |20| . This difference can be explained 
by a 7% offset in the energy calibration alone, well within 
the stated uncertainty of the two experiments. 

The GZK effect predicts a suppression in the UHECR 
flux above 10 19 ' 8 eV. We fit our data to a model consist- 
ing of galactic and extra-galactic sources |2l|, that in- 
cludes the GZK effect. We use the extra-galactic source 
model of Berezinsky et al. [2^] , where we assume a uni- 
form (over the Universe) proton source distribution with 
a maximum at-source energy of 10 21 eV, and a galactic 
spectrum consistent with observations that the composi- 
tion changes from heavy to light near 10 18 eV. The x 2 
of this fit is 48.5 for 37 degrees of freedom, and the fit 
is shown in Fig.0] Details can be found in jilij . In this 
model the fall-off above logi? of 19.8 is due to cross- 
ing the pion production threshold, and the second knee 
comes from e + e~ pair production pile-up. 

For comparison, the published AGASA spectrum 0] 
is shown in Fig. Compared to the HiRes monoc- 
ular spectrum, the AGASA flux values are about 60- 
70% higher. The AGASA data contain 11 events above 



10 20 eV and 24 above 10 19 8 eV, with an integrated ex- 
posure of ~ 1.0 x 10 3 km 2 sr-yr. The HiRes data have 
two events at or above 10 20 eV, and 10 above 10 19 8 eV, 
with exposures of 2.4 and 2.2 x 10 3 km 2 sr-yr. at these 
two energies. If we were to increase the HiRes energy 
scale by 17% (one RMS systematic deviation), the num- 
ber of events above 10 20 , and 10 19 8 eV would become 
three and 20. A decrease in energy scale of 17% changes 
these numbers to one and 6. 
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